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We propose an efficient protocol for braiding atomic Majorana fermions in wire networks with 
AMO techniques and demonstrate its robustness against experimentally relevant errors. Based on 
this protocol we provide a topologically protected implementation of the Deutsch-Josza algorithm. 



The prediction of particles with anyonic statistics in 
topological phases of matter has resulted in the proposal 
of decoherence-free Topological Quantum Computation 
(TQC) [T}]3]. TQC requires the creation of anyonic par- 
ticles as well as their controlled interchange, known as 
braiding, which is the fundamental building block of topo- 
logical quantum gates ^ [5] . While the implementation 
of these tasks in real physical systems is an outstanding 
challenge, the reported observation of anyonic Majorana 
fermions (MFs) in hybrid superconductor-semiconductor 
nanowire devices [BHH] and the proposals for the manipu- 
lation [9l-[TT] of anyonic Majorana fermions (MFs) in solid 
state systems are promising first steps in this direction 
|llH16j . A complementary and promising approach to- 
wards realizing and coherently control MFs are ultracold 
atoms confined to ID optical lattices coupled to BCS or 
molecular atomic reservoirs. The recent realization of a 
quantum gas microscope [171 HH] for optical lattices adds 
single site addressing and measurement to the toolbox 
of possible atomic operations to create and detect MFs 

[TMH- 

Building on these experimental advances, we describe 
in this Letter an efficient braiding protocol for atomic 
MFs, based on performing simple lattice operations on a 
few sites in an array of ID wires, and we provide a careful 
study of the full braiding dynamics including imperfec- 
tions. In addition, we will show that these elementary 
braiding operations, although they do not represent the 
complete set of quantum gates [551 123] , can be combined 
to realize a Deutsch- Jozsa algorithm [21] , demonstrating 
that the implementation of simple quantum algorithms 
in atomic topological setups is within experimental reach. 

Braiding of atomic Majorana fermions. We consider 
a system of single component fermions that are con- 
fined to an array of one-dimensional (ID) wires of L 
sites (see Fig. [ij and that are governed by a Hamil- 
tonian H = En^^"-*- The Hamiltonian = 

-J^n^janj+i + Aan,jan,j+i + h.c. - fiJ2n ^ti^n re- 
alizes a Kitaev chain [5S] in the n-th wire. The operators 

a|j J and Onj are fermionic creation and annihilation op- 
erators, J > and A G E, are nearest-neighbor hopping 
and pairing terms, and is a chemical potential. As 
demonstrated in [19], a Hamiltonian of the form i?'"-' 
allows for a cold atom implementation: While the hop- 




FIG. 1: Realization of an array of one-dimensional Kitaev 
wires in an optical lattice setup: Atoms (red circles) can hop 
between neighboring sites (blue circles) with strength J along 
the individual wires. The pairing term of strength A can be 
realized by a Raman induced dissociation of Cooper pairs (or 
Feshbach molecules) forming an atomic BCS reservoir. 



ping term arises naturally in an optical lattice setup, the 
pairing term can be realized by a Raman induced disso- 
ciation of Cooper pairs (or Feshbach molecules) forming 
an atomic BCS reservoir. 

It has been shown in [5S] that the Hamiltonian i?*^") 
supports zero energy Majorana fermions of the form 
^''l/r. ~ j^c„ J- with (real) coefficients v^^f^ which 

are localized at the left/right end of the n-th wire. Here, 

Cn,2j-i = alj + On J and Cn,2j = (-i)(alj - "nj) are 
Majorana operators fulfilling {cn.k, Cm.i} = 2SkiSmn- For 
the "ideal" quantum wire (J = |A|,/i = 0), one has 

v^i — '^,v^2L — 1 ^^'^ ^1^^ ''^n^f ~ Otherwise, 
the modes 7^"'*^ decay exponentially inside the bulk. 
Each wire has two degenerate ground states |0„) and 
|1„) with even and odd parity, respectively, correspond- 
ing to the presence or absence of the Majorana fermion 
fn = Jl '' - ilR\ i-e. /„|0„) = 0, /t|0„) = |1„). For a 
proposal how to prepare MFs in the desired parity sub- 
space see [20] . 

Since MFs exhibit anyonic statistics, an appropriate 
interchange of two Majorana modes, 71,2 1 allows to real- 
ize the braiding unitary Ui, = £^^7172/4^ This interchange 
which leads to the transformation 71 i-> —72, 72 H> 71 
resulting in a non-trivial phase factor for the wave func- 
tion is the key step for realizing a TQC. In the following 
we present a protocol for a cold atom implementation 
that allows to realize braiding. To this end, we con- 
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FIG. 2; Braiding protocol for two perfect quantum wires. The 
zero-energy Majorana modes that are initially on the upper 
(lower) wire are shown as yellow (black) spheres, while the 
blue ones corresponds to the Majorana operators which are 
coupled into finite-energy fermionic modes. Coupling of Ma- 
jorana operators via hopping and pairing (Kitaev coupling) is 
indicated by grey solid links, while the coupling via hopping 
only is shown as a dashed link. 



sider two neighboring wires n and n + 1 governed by two 
ideal Kitaev Hamiltonians and The use of 

ideal wires allows for a simple analytic treatment because 
only six Majorana operators are involved in the proto- 
col, ft is convenient to label the sites by {w,j), where 
w = u,l denotes the upper (n) resp. lower (n + 1) wire 
and j = 1,. . . ,L. We label the sites that are involved 
in our protocol as si = (w, I), S2 = (u, 2), S3 = {1, 1) and 
S4 — {1,2) (see Fig. [2]). To simplify notation, we write 

Cu.j = Cj, and cij = dj with Majorana modes j^^"^ = Ci, 

Let us now show how to braid the left Majorana modes 
7^"'' and around each other with only local (adia- 
batic) changes in the Hamiltonian on the left edge of the 
system. These changes include switching on/off (i) the 
hopping i?^''^^-. = —Jai ag^ + h.c. and (ii) the pairing 

Hi^K — Jag^ag- + h.c. between the neighboring sites Si 

Alp) 



between sites Si — S2 and S3 — S4, and, at the same time, 
switch on the hopping between sites si — S3: 

Hi{i) = cosMH^^- + ) + sin (l,tHi''l 

-* V / S1S2 ^ S1S3 

= -iJ [cos0t(c2C3 -I- ^2^3) + sin(/>t(c2di - cic?2)/2] . 

During this process the zero modes evolve accord- 
ing to j'i\4>t) = {2cos(f>tCi - sin(/)td3)/v/l -I- Scos^ ^j, 
7L^(0t) = (2cos0tdi - sin(/)tC3)/yr+l5cos2^, such 
that at the end 7^^"'' = — ^3 and j^^^ = — C3. Note that the 
two decoupled sites si and S3 carry exactly one fermion 
which has been taken out of the system. 
Step II: We put now this fermion in the lower wire by 
switching on Hi^l^ between sites S3 — S4, and H^\^ be- 
tween the sites si — S3: 



Hii{t) - i/lH + sin ( H^^l + i/l^J 

-'-'V / S1S3 T-c S1S3 S3S4^ 

- if [{c2di ~ £1^2) + sin<j)t{c2di + Cid2 + 2^2^3)1 • 

The zero modes evolve as J^^\4>t) — (2sin0(Ci — (I — 
sin(/>t)d3)/-^4sin^ 0t + (I - sin(/)t)2, = -C3, such 

that at the end 7^^"^ = ci and ^'"^'^ — — C3. Note, that 

at this stage the Majorana mode 7^^"'' (7^^''') has already 
been moved from the upper (lower) to the lower (upper) 
wire. However, two additional steps are needed to re- 
cover the original configuration of the wires. 
Step III: We move the Majorana mode from the site s\ 
to the site S3 by switching on H^g"^^ and simultaneously 
switching off -ff^^'' between the sites si — S3 : 

HuM = sin0,i/f ) -f cos -f 42 

= —iJ(c2d\ cos (fit + (^2^3) ~ iV^sin0tCiC2. 



The evolution of the zero mode j^j^^ = (Jcos^tCi -I- 

V sin (j)tdi) / {J cos cft)"^ + (ysint/it)^ results in 7^^"' — 
(I) 



Step IV: Finally, we switch off Hi^^^ and switch on Hi^J^ : 

Hjy it) = Sin ^tHifl + ilg + cos cl>tHi[^^ 

= —iJ [sin 4>tC2C3 + (i2C?3] — iV cos 4>tCiC2 



The zero modes are given by 7)^ 



and Sj, and (iii) the local potential Hi^' = 2Valag. on di, while 7^'^ = -C3 remains fixed 
site Si. Note that a combination of (i) and (ii) allows to 
switch on/off the Kitaev coupHng Hi^l, = Hi'^\, . 
These operations are based on the single site/link ad- 
dressing available in cold atom experiments [17l [26] . 

Let us now describe the braiding protocol in detail. 
The physical process behind is the transfer of one fermion 
from the system (i. e. either from the upper or from the 
lower wire) into the lower wire. We characterize the re- 
quired adiabatic changes via a time-dependent parameter 
(j)t that varies from to 7r/2, and perform them in four 
steps. In describing these steps, we will only write down 
the Hamiltonian for the four involved sites and follow the 
evolution of the zero modes which are always separated 
by a finite gap from the rest of the spectrum. 

Step I: We decouple the two very left sites, si and 
S3 from the system by switching off the couplings Hi^l 



in) 



.1 (0 



so 



-(Jsin^tci -I- ycos^tC3)/v'(Jsin0t)2 -|- {V cos (pt)^ 
that finally we get the desired braiding 7^"' i~> jf' and 

j^l^ 1-^ ~7l''' of left the Majorana modes on the wires n 
and n -|- 1, which corresponds (up to unimportant phase 
factor) to the unitary J7„ — e^'^L ''''i'/'*. 

Note that the braiding in the other direction, Ul and 

7^^"'' 7^''' I— >■ 7i'\ can be achieved by putting 

the uncoupled fermion in the upper (instead of the lower) 
wire with a simple modification of Steps II-IV. 



3 




FIG. 3: Evolution of the Majorana correlation functions 
(ni''7i'|) (red, o), {z7i''7l^") (blue, v), (ni'^T^,") (red, □), 
and {iTi^'lj^') (blue. A) during the braiding protocol with 
errors a in the local operations for two non-ideal quantum 
wires with |A| = 1.5 J and /x = 0. Markers are only drawn in 
regions where the correlation functions are non-zero. 
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The braiding results in the change of the correlation 
functions of the Majorana operators (see Fig. 2) and thus 
changes also the long-range fermionic correlations. This 
can also be translated into the change of the fermionic 
parities of the wires: If |0„) (|ln)) denotes the state of 
the n-th wire with even (odd) parity and, for example, 
we start from the state |0„0„+i) with both wires with 
even parity, then the braiding [/„ results in C/„|0„0,i+i) — 
) + |l„l„+i))/\/2, and C/2|0„0„+i) = |l„l„+i). 
The result of the braiding, therefore, can be checked by 
measuring the change of the Majorana correlation func- 
tions in Time-of-Flight or spectroscopic experiments |20j , 
or by measuring the parity of the wires by counting the 
number of fermions modulo two !l8j. 

Non-ideal wires and non-perfect operations. We have 
just demonstrated the braiding for the case of ideal Ki- 
taev wires and perfect local operations (single site/link 
addressing). Remarkably, the topological origin of the 
Majorana modes ensures the robustness of the results of 
the braiding protocol based on Steps I- IV also in the real- 
istic case of non-ideal wires and local operations provided 
the Majorana modes are spatially well-separated. We 
have checked this numerically by considering two non- 
ideal wires with J ^ |A|, ^ 7^ and assuming that the 
local operations have an error a in the following sense: 
(i) Switching on the hopping J and/or the pairing A 
between the sites [u, 1) — (d, 1), also introduces the hop- 
ping Ja and/or the pairing aA between the adjacent 
sites (m, 2) — (d, 2). (ii) Switching off the couplings be- 
tween the sites (w, 1) — {w, 2) also reduces the couplings 
between the sites (w, 2) — (w, 3) by a factor (1 — a), (iii) 
Raising the local potential V on the site (u, 1) results in 
a local potential aV on the neighboring sites (u, 2) and 
{1, 1). As an example, we present in Fig. [s] numerical re- 
sults of the braiding protocol with errors a — 0.05 and 
a — 0.1 in the local operations for two quantum wires 
of the length L = 40 with |A| = 1.5J and /x = 0. One 
can clearly see the robustness of the final results of the 



FIG. 4: Braid group in a setup of three wires. We present the 
real-time evolution of the correlations functions *(7i"''7/{"') 
under the action of {a)U2U\, (b) U1U2 and (c) U1U2U1 for a 
chain of the length L = 40 with | A| = 1.5J and fj, = 0. Mark- 
ers are only drawn in regions where the correlation functions 
are non-zero. 



braiding. 

Braid group. It is also easy to check that the uni- 
tary transformations C/„ of the Majorana operators cor- 
responding to the braiding protocol fuUfill all necessary 
conditions of the braid group [2] : For any two braiding 
unitaries U„ and Un+i, one has C/„C/„+i 7^ C/„+iC/„ and 

Un-lUnUn-l = UnU„-iUn. 

To demonstrate this, consider three wires with left Ma- 
jorana modes 7^^^"*, 7^^'* and 7^^'*, and braiding unitaries 
Ui = e^^L^iTl^ and = e'^^i'Ti'/^ that braid the 
modes 7^^^'*, 7^^'' and 7^ , \ respectively. The braid 
group conditions for U\ and Ui then immediately follow 
from the following formulae 



jj jj I (1) (2) (3)n 



(2) (3)x 

7l ,7l ) 



(7f, 
(7i^\ 



-7i^\ 



-Tl ) 



(2) (1)n 

-7l ,7l ) 



(1) 



(2) 



These properties can be tested experimentally by mea- 
suring the corresponding changes of the fermionic corre- 
lation functions. For example, the action of U1U2 results 

• •/ (1) (1)\ . ■/ (3) ■/ (2) (2)\ ^ ., (1) (2)> 

and iil'^^'l^R) ^ ^''■{"/l^Ir^)^ while U1U2U1 pro- 
duces the following changes: ii'^^^jR^) '-^ '''{'Il^Ir^)^ 

■I (2) (2)\ , ■; (2) , ., (3) (3), ^ ., (1) (3), 

ATl Tr ) ^ -^ilL Tr ), and i{lh Tr ) ^ ^ilh Tr ) 
(see Fig. [4]). This change in the correlation functions 
can be measured, for example, in TOF or spectroscopic 
experiments as proposed in Ref. (20j . 



4 



Deutsch-Josza algorithm. Although the braiding of 
MFs is robust, it does not provide a tool to construct 
a universal set of gates needed for TQC: As it has been 
shown in Ref . P3j , only a subgroup of the Clifford group 
can be realized via braiding. Fortunately, not all QC al- 
gorithms require a universal set of gates. One example is 
the Deutsch-Josza algorithm [53] which, as we will show 
below, can be implemented for two qubits in a remark- 
ably efficient way via braiding of MFs. 

The Deutsch-Josza algorithm allows to determine 
whether the function ("oracle") g{x) which is defined on 
the space of states of n qubits and takes the values or 1, 
g : {|0) , I— {0, 1}, is constant (has the same value, 

say, 0, for all inputs) or balanced (takes value for half of 
the inputs, and 1 for the other half). For the algorithm 
to work, the function g has to be implemented as the 
unitary Ug : ^ (-l)f(^)|x), where \x) G {|0) , 
which is actually a major problem for experimental re- 
alizations; A faulty oracle spoils the quantum speedup 

For two qubits with the computational basis 
{|00), |01), |10), 111)}, a possible choice for Ug is 



g, - diag(l, 1, 1, 1), Ug, = diag(l, 1,-1, -1) 



diag(l, -1,-1,1), Ug^ =diag(l, -1,1,-1), 



for the constant go and the balanced gi, 52 , and 173 oracle 
functions, respectively. (Note that an equivalent set of or- 
acles can be obtained by multiplying the above unitaries 
with —1.) The algorithm works then in the following way: 
After preparing the system in the state |00), we apply the 
Hadamard gate H to each qubit, H \0) = {\0) + \1))/V2, 
H\l) = (|0) - \l))/V2, then we apply the unitary Ug 
corresponding to the oracle under test, then again the 
Hadamard gate to each qubit, and, finally, we measure 
the probability to find the system in the state |00). This 
probability is 1 if g{x) is constant, and if g{x) is bal- 
anced, as can be seen from the following calculations 

ioo)"^"^5:ix)j^^^(-i)^(>')ix) 

X X 



H^H 1 



where we define x — {xi,X2) and y = (2/1,2/2) with Xi, 

€ {0, 1}, and X • y = xiyi + X2?/2- 

To implement the above algorithm, we use a setup of 
three quantum wires in the geometry shown in Fig. [5] 
and define a computational basis for two qubits as 



100) 



A\Of), |01) = /]|0/), |10) = fl\Of), and 



|11> = /1V2VIIO/) [23- Here |0/) is the vacuum state 
for fermionic modes fi = (7!^'' — i^[^^)/2, where ^[^^ 
and 7j are two Majorana modes on the i-th wire. 
Note that in this setup with three wires we encode only 
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FIG. 5: a) Setup and implementation of the "oracle" Deutsch- 
Josza algorithm for two qubits via braiding, b) Implementa- 
tion of the "oracle" unitary Ug via braiding (see main text). 



two qubits [28 . This is because the braiding preserves 
fermionic parity and, therefore, all states from the com- 
putational basis must have the same parity (odd in our 
case). 

The Hadamard gates and the oracle unitaries Ug^ can 
be implemented by noting that the braiding of Ma- 
jorana modes 7,; and 7j is equivalent to the unitary 



u 



(1 



7i7j)/-\/2. Then, it follows imme- 
diately that H®t = U12U23U12 and t(E)H = U^^U^^U^^i 
for the Hadamard gates acting on the first and the sec- 



ond qubit, respectively, and Ug-, = U(2^ U, 



•^56 



92 



- TT"^ 

— "^34; 



and 



for the oracle unitaries {Ug^ — 1). As a re- 
sult, the Deutsch-Josza-algorithm can be realized with 14 
braiding operations. In our case, however, the number of 
operations can be reduced to nine: The sequence 



= Ui5U5eU23Ui2U,^U56U45Ui2U: 



^23 



acting on |00) gives Z^o|00) — |00) for the constant case 
and Ui\00) = i\10), Z^2|00) = |11), and Z^slOO) = i|01) 
for the balanced case. Note also that this protocol can 
be implemented in five steps because operations on the 
Majorana modes 71,2,3 and 74.5,6 before and after the or- 
acle unitary Ug. can be performed in parallel. The final 
state of the system and, therefore, the probability to find 
it in the state |00), can be determined by measuring the 
parities of the individual wires in a spectroscopic experi- 
ment |20| or fermionic number counting |18j . Taking into 
account the discussed insensitivity of the braiding to ex- 
perimental imperfections, the proposed protocol provides 
a robust implementation of the Deutsch-Josza algorithm. 

Conclusion. We have presented an efficient way of 
braiding MFs in a cold-atom setup and used it to im- 
plement the Deutsch-Josza algorithm in a topologically 
protected way. By adding well-controlled though topo- 
logically unprotected operations (e.g. the SWAP-gate), 
one can go beyond the braid group and provide a univer- 
sal "hybrid" set of gates for quantum computation (see 
also Ref. J29l |30j). We address this issue in our future 
work. 
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